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Abstract 

Pauli exclusion between the carriers of ./V excitons induces novel many-body 
effects, quite different from the ones generated by Coulomb interaction. Using our 
commutation technique for interacting close-to-boson particles, we here calculate 
the hamiltonian expectation value in the iV-ground-state-exciton state. Coulomb 
interaction enters this quantity at first order only by construction ; nevertheless, due 
to Pauli exclusion, subtle many-body effects take place, which give rise to terms in 
(Na^,/V) n with n > 2. An exact procedure to get these density dependent terms is 
given. 
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It is known that excitons being made of two fermions are not exact bosons. Their 
underlying fermionic character is however a major difficulty which has been overcomed 
very approximately up to quite recently. An approach C 1 ' 2 ) which has looked reasonable 
for years is to assume that excitons are exact bosons provided that their close-to-boson 
character is included in an appropriate exciton-exciton interaction which is basically a 
Coulomb interaction dressed by exchange processes. 

Very recently, we have developed a formalism which allows to include the close-to- 
boson character of the excitons exactly, through an extremely simple and physically quite 
transparent algebra ®. Using this "commutation technique", we have already shown 
( 4 ), by calculating the correlations of two excitons, that Pauli exclusion between the two 
electrons and the two holes of these two excitons, enters their Coulomb terms in such a 
subtle way that a nai ve bosonic hamiltonian for excitons t 1 ' 2,5 ) cannot produce these terms 
correctly beyond first order in Coulomb interaction, reducing considerably the impact of 
such an effective hamiltonian. 

A way to physically understand the problem is to realize that excitons feel each other 
through both Coulomb interaction and Pauli exclusion between their components. Con- 
sequently, besides the usual many-body effects resulting from Coulomb interaction, the 
excitons do have very unusual ones coming from Pauli exclusion. These two kinds of 
many-body effects being basically independent, they enter the various quantities of phys- 
ical interest quite differently. So that there is no reason, for the exchange processes 
resulting from Pauli exclusion, to dress the Coulomb interaction in an unique way. In 
other words, there is no reason to describe the physics of interacting excitons by one 
dressed exciton-exciton interaction only, as assumed in the effective bosonic hamiltonian. 

In the present work, we calculate the semiconductor hamiltonian expectation value in 
the iV-ground-state-exciton state. We here forget about the spin degrees of freedom ( 6 ) of 
the excitons, for the sake of simplicity. In the absence of Coulomb interaction and Pauli 
exclusion, this hamiltonian expectation value should be equal to NE with E Q being the 
ground state energy of one exciton. Due to Coulomb interaction, which by construction 
enters this quantity at first order only, we should find one additional term resulting from 
the average Coulomb energy in this iV-exciton state, and which should be Na x /V smaller 
than the main energy NE Q , a x being the exciton Bohr radius and V the sample volume. 
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The exact calculation shows that the hamiltonian expectation value also contains terms 
in (Na>x/V) n with n2, which originate from many-body effects between the N excitons 
induced by Pauli exclusion. Using our commutation technique, we show how we can derive 
this iV-exciton average energy at any order in r\ = Na^/V in a systematic way. 

The paper is organized as follows. In section 1, we recall the basic relations of our 
commutation technique and use them to calculate the expectation value of the hamiltonian 
in the iV-ground-state-exciton state, < H >jv. In section 2, we derive the recursion 
relation satisfied by the matrix elements which appear in the expression oi < H > N , and 
we iterate this relation to obtain the density leading terms of these matrix elements. In 
section 3, we calculate the three first terms of the 77 = Nal/V expansion of < H > N 
explicitly. In the conclusion, we discuss the obtained result and relate it to the new 
criterion for bosonic behavior of excitons we recently proposed ( - 7 \ 

1 Expectation value of the hamiltonian using our 
commutation technique 

The semiconductor hamiltonian expectation value in the iV-ground-state-exciton state 
reads 

< H >K= < wg(gi)> > t (1) 

< t>|(B„HBj» > 

where |i> > is the electron-hole pair vacuum state and Bq the (exact) creation operator 
of one ground state exciton, i. e. exciton with a center of mass momentum Qo = and 
a relative motion in its ground state ip UQ . The algebraic calculation of < H > jv is quite 
easy through our commutation technique. 

Two important relations of this commutation technique are ( 3 ) 

H, Bj] = E t Bj + V} , (2) 

V},B^=Y.C nij BlBl , (3) 

ran 

where E { is the % exciton energy and B\ the (exact) creation operator of an % exciton, % 
standing for [y^ Qj). The operator V} comes from the Coulomb interactions between the 
i exciton and the rest of the system, while the parameter ^ r ni j corresponds to the direct 
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Coulomb scattering of the excitons into the (m,n) states. In r space, ^nij reads ( 3 ) 

Zmnij = £jmn = \ J dr e, dv hl dr e2 dr h2 (\>*Jt^ ^l)0n( e 2, h 2 ) 

x [V ei e 2 + V hlh2 - V eih2 - V e2hl ] 0i(ei, h 1 )(j) j (e 2 , h 2 ) + (i <-> j) , (4) 

where <f>i(e, h) is the total wave function of the i exciton, i. e. the product of the relative mo- 
tion wave function (p Vi (r e — r^), and the center of mass wave function e l Qi-( a er e +a h r h ) j 
with a e = 1 — ah = m e / (m e + rrih)- 

Using Eqs (2,3), it is straightforward to prove by induction that 



H, (Blf 



so that, as H\v >= = Vq\v >, we get with H acting on the left 

NF | N(N-l) dir < v\{B ) N - 2 B m B n (Bl) N \v > 

1 mn < V\(B ) N (B^y y \V > 

As expected, Coulomb interaction enters < H > N at first order only, as < H > N is linear 
in the Coulomb scatterings £oomn- 

If the excitons were exact bosons, the ratios in the sum would be 1 for m = n = 0, 
and otherwise, so that we would get < H >n= NEq + N(N — 1)£qqq /2. We can then 
note that for n — % — j, Eq (4) reduces to 

& = y 2 j ^ dr hl dr e2 dr h2 ^-Q-M^+a^)^^ _ rfc j ^(^ _ r ,J 

N|2 / e 2 e 2 e 2 e 2 \ 
x \ l Pu l { r e 2 - r h2 )\ I - + - --- --- 1 

\|*"ei 1 \^-hi Yfo 2 1 |r ei r^ 2 1 |I" e2 Y^] J 

= 0, (7) 

as can be seen by exchanging r e2 and r^ 2 . Consequently, strangely enough, if we forgot 
about the close-to-boson character of the excitons, < H > N would be equal to NE as if 
the excitons were not interacting at all ! 

Of course, excitons are not exact bosons. Let us now show how we can calculate 
the matrix elements appearing in the sum of Eq (6) for exact excitons. Here again, 
two important relations of our commutation technique make their calculation quite easy, 
namely ® 

' n = %-A, , (8) 
2£A mmjJ Bt. (9) 



B h B) 



D B 1 

"mil ■ LJ j 



Dij is the deviation-from-boson operator, while the parameter X m nij describes the fact 
that there are two ways to couple two electrons (ei,e 2 ) and two holes (hi,h 2 ) to make 
two excitons, either (ei, hi)(e 2 , h 2 ) or (ei, h 2 )(e 2 , hi). In r space \ m nij reads ® 

Ammj = Kjmn = \j dr e 1 dr hl dr e2 dr h2 (f)* m (e 1 , h)(f)* n (e 2 , /i 2 )0i(ei, h 2 )<l) j (e 2 , h x ) 

(10) 

Eqs (8-10) are the key equations which allow an easy calculation of the matrix elements 
appearing in < H 

2 Calculation of the matrix elements appearing in 

< H > N 

From Eqs (8,9), it is straightforward to prove by induction that 

B m , (Bl) N ] = NSmoiBl)"- 1 - N(N - 1) ]T A^flJ^-a _ N (Bl) N ^D m0 . (11) 

v 

Before calculating the matrix elements of < H >jv, with two excitons possibly outside 
the ground state, let us first consider the simpler ones, with only one exciton outside the 
ground state, namely 

A^ =< „ | {B,) N - l B m {Bl) N | v > /N\ . (12) 
As | v >= 0, we get from Eq (11) that A^ verifies the recursion relation 

= 5 mQ A[ N ~ l) - (iV - 1) E A mP oo {Af-^y . (13) 

p 

If we iterate this equation, and note that Aq^ is nothing but the quantity F N defined in 
Ref. (7), 

A ( N) = F N =<v | (B ) N (Bl) N | v > /N\ , (14) 
we find that A^P expands as 

= F N -i8 m0 — (N — l)F N - 2 \ m000 + (N — 1)(N — 2)F N _ 3 E A mp ooAoopo 

p 

_( iV _i)(iV_2)(iV-3)F J v_4E A ™pooAoo P Aooo + -- - • (15) 

pq 
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By expliciting the A's, it is possible to show that, for m — 0, Eq (15) is nothing but the 
recursion relation between the i*Vs given in Ref. (7), namely 



F N = F N ^ - (N - 1) a 2 F N _ 2 + (N-1)(N- 2) a 3 F N _ 3 , (16) 



with 



°n = El $ -o( k ) 



2n 



in 3D, $^ (k), given in the appendix, being the Fourier transform of ip UQ (r). 

Let us now turn to the matrix elements with two excitons possibly outside the ground 
state, namely 

B£2 =< v | (B ) N - 2 B m B n (Bl) N I v>/N\. (18) 

From Eq (9), we get 

D mi , (Bl) N ] = 2Nj2^oBl(Bl) N - 1 , (19) 

q 



easy to prove by induction. Using Eqs (11,18,19), we find that B^J verifies the recursion 
relation, 

&mn = F N _ 2 {^mQ^nQ — A mri oo) 

- (N-2)J2 (A ( q N ~ 2) ) * | ^m0A„ 9 00 - E X mp00K q p^j + ipi <-> n) j 

+ (N - 2) (TV - 3) ]T (C" 2) )* A mp00 A ng00 . (20) 

P<3 

We can then use Eq (15) for *4.^ v - ) and iterate Eq (20). The three first terms of the 
expansion of B^J are obtained by replacing by the two first terms of Eq (15) and 

B iq ~ 2) b y the first term of E q (20). This leads to 

Bmn = ^V-2(^m0^n0 ~~ A mn0 o) 

-(N - 2)Ftv_ 3 I ^ m0 A„ooo - X! ^mpooAnOpo^j + (m <-> n) j 

<5m0 X/ -^npOoAoOpO ~~ X] ^mpOO AnqpoAoOgO J + ( m n ) 
\ P P<J / 

+ A m oOoA n OOO — X/ ^mpOoAngOoAoOpg \ + " ' ' (21) 
P<3 J 



+ (iV-2)(iV-3)F w _4 



The first term of B^J corresponds to all Pauli couplings between two ground state excitons 
(0, 0) and the two (to, n) excitons; the second term of B^ r ] corresponds to all Pauli 
couplings between three ground state excitons (0,0,0) and the three (to, n, 0) excitons; 
the third term couples the (0, 0, 0, 0) excitons to the (to, n, 0, 0) excitons, and so on . . . 

3 Density expansion of < H > N 

Let us now return to the expectation value of the hamiltonian < H Using Eqs 

(6,14,18,21), and the fact that £ooom =0 ( see E Q ( 7 ))> we find that < H > N reads 

< H > N = N(E + A) , (22) 

where the three first terms of A are given by 

- iN -l)( N -2)(N-3)*j^(s s + % -§)+"• , (23) 

the expressions of the various sums S n in terms of £ dir 's and A's being given below. 
By using Eqs (4,10), it is easy to see that the first sum Si reads in r space 

mn 

= J dr ei dr hl dr e2 dr h2 <f>o(e u ^i)0S( e 2, he) 

x (V eie2 + V hlh2 - V eih2 - V e2hl )(j)o(e 1: h 2 )(f) ( e 2, h) , (24) 

which is nothing but £oooo = ^oooo defined in our previous works dealing with the matrix 
elements of H in the two-exciton subspace This term is also the usual exchange 
Coulomb term of the effective bosonic hamiltonian for excitons ( - 2 ' 8 \ 
The second sum of Eq (23) involves three excitons. It reads 
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^2 — ^ £oOmn^mpO(AnOpO 
mnp 

= J dr ei dr hl dr e2 dr h2 dr e . A dr h . A 0*(d, h)(j)* (e 2 , h)(j)* (e 3 , h 3 ) 

x [V eie2 + V hlh2 - V eih2 - V e2hl ] o (ei, h 2 )(f) (e 2 , h 3 )0o(e 3 , h) • (25) 

The three sums of the third term of Eq (23) involves four excitons. The first sum S3 
reads 

S3 = ^2 ^OOmn^rnpOoKtqpO^OOqO 
mnpq 

= J dr ei dr hl dr e2 dr h2 dr e . A dr hs dr e4 dr hi 0* (d , hi ) (j)* (e 2 , h 2 ) 
x 0o( e 3, h 3 )(f)* (e 4 , h 4 )[V eie2 + V hlh2 - V eifl2 - V e2hl ] 
X(j>o(ei,h 2 )<f>o(e2,h 3 )(f>o(e3,h4)(f>o(e4,hi) . (26) 

The sum S' 3 , defined as S' 3 = T,mn P q Coomn^mpOoKqOO^oopq , reads as Eq (26) except for 
the four O ' S which are replaced by 

0o(ei, h 3 )4> (e 2 , h 4 )<i) (e 3 , h 2 )<f) (e 4 , h) . 

The last sum S 3 , defined as S'3 = J2 mn £oomn^moooA n ooo, reads as Eq (26) except for 
the four o 's which are now replaced by 

0o(ei, h 3 )4> (e 2 , h 4 )<j) (e 3 , ^i)0o(e 4 , h 2 ) . 

This sum is in fact equal to zero as easily seen by exchanging (r e2 <-> r h2 ), and (r e4 <-> r^J, 
in the integral. 

It is physically important to note that all these sums contain the same Coulomb 
coupling between only two excitons made with (ei,hi) and (e 2 ,h 2 ). Couplings between 
more than two excitons appearing in S n2 , result from many-body effects induced by the 
close-to-boson character of the excitons. 

The sums Si, S 2 , S 3 , S' 3 are calculated in the Appendix (see Eqs (40-43)). Although 
not obvious at first, Si, S 2 and S 3 + S 3 /2 are in fact real even if $> vo is not, as necessary 
to have a real energy change A induced by these Coulomb and Pauli couplings. 

A also contains density dependent terms through the ratios Fn- p /Fn which differ 
from 1 due to the close-to-boson character of the excitons. These ratios can be obtained 
from the recursion relation (16), which also reads 

-p- = l-(N-l)a 2 ^ + (N-l){N-2)a 3 ^ + -- (27) 

" N—l r N _i -Fn-1 
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As a n is in (a? x /V) n 1 (see Eq (17)), we can iterate Eq (27) to obtain the ratio Fn/Fn-i, 
for TV >> 1, as an expansion in powers of the parameter 77 = Na^/V. We get 

Jk_ = 1 _ Na2 + iv2 ((T3 _ a 2) + {rf) (28) 

F N—l 

As for large TV, Fn- p /Fn — {F^-i/ Fn) p , the ratios appearing in Eq (23) are given by 
^ = l+2Na 2+ N 2 (5al-2a 3 )+0(r ] 3), = l+3Na 2 +Otf) , *£=i = l+Ofa) . 

-TAT -TAT ^JV 

(29) 

By inserting Eq (29) into Eq (23), we obtain the following expansion of the energy change 
A in powers of 77 : 



A = -N^- + N 2 [S 2 - (72^] + N 3 



-S 3 -^ + 3<j 2 S 2 -^(5<j 2 2 -2<j 3 ) 



+ 0( V 4 ) (30) 



Using Eqs (17,40-43) and putting everything together, we finally find that the hamilto- 
nian expectation value in the iV-ground-state-exciton state expands in powers of 
i] = Na 3 x /V as 

im / 13tt 73tt 2 2 3517tt 3 3 4 .\ 
<H> N = NE (l-— V * + 0(riy 

= NE (l - 13.61 7^ + 36.02 if -519.3 rf + 0(rf)) (31) 

Let us again stress that terms in rf 1 with n > 2 come from many-body effects induced by 
Pauli exclusion, since Coulomb interaction enters < H > N at first order only. 



4 Conclusion 

We have calculated the expectation value < H > N of the exact semiconductor hamilto- 
nian in the iV-ground-state-exciton state (Bq) n \v >, using the commutation technique 
we recently introduced. Due to novel many-body effects induced by the close-to-boson 
character of the excitons, < H >n appears as an expansion in powers of the density 
through 77 = Na^/V, Coulomb interaction entering this quantity at first order only by 
construction. Higher order terms in rf 1 , with n > 2, result from both, sophisticated ex- 
change processes in which the Coulomb interaction appears at first order only, and purely 
Pauli many-body effects which make < v \(B ) N (Bq) n \v > to differ from its exact boson 
value A!. 
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From the result given in Eq (31), we see that the prefactors of the expansion of < H >n 
are rather large so that r\ = Na^/V has to be much smaller than 1 for < H >n to be 
equal to the energy of N non-interacting boson-excitons NE . In a previous work ( 7 \ 
we have already shown that, while the Mott criterion corresponds to rj = Na^/V ~ 1 
for the electron-hole pairs to be bound in excitons, i. e. for the excitons to exist, the 
criterion for the excitons to behave as bosons is more like lOONa^/V ~ 1. As the 77 
expansion of < H > N is physically linked to many-body effects induced by the close-to- 
boson character of the excitons, it is after all reasonable to find similar conditions for 
excitons to behave as bosons and for the hamiltonian expectation value to be equal to the 
one of N non-interacting boson-excitons. 
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APPENDIX 



We start from the expressions of Si, S2, S3 and S3 in r space given in Eqs (24-26), 
and we rewrite them by using the Fourier transforms of the exciton wave function and 
Coulomb potential. For Q = 0, we have 



e 
r 



(32) 
(33) 



with Vq = 4vre7Vg 2 in 3D. 

If we insert Eqs (32,33) into Si given in Eq (24), and we perform the integrals over 
the r's, we get 



l<Mk' 



■'- 2 -<t>: (k)<uko 

By using the Schrdinger equation for the exciton relative motion, 



Si = 2^y k _ k # J , (k) 

k,k' 



(34) 



(h 2 k 2 /2f, - e^)^(k) - £ V *-v $ -o(k') = , 



(35) 



with e uo — E — —e 2 /2a x , we can check that Si is real. 

If we now insert Eqs (32,33) into S2 given in Eq (25), and perform the integrals over 
the r's, we get 



s 2 = e fk-* f2|$, (k)i 4 i$, (koi 2 - i$, (k)i 4 $: (k) $ vo (ko 



k,k' 



- «(k)<f, () (k)) ($ 2 (k')$: (k' 



(36) 



The second term of Eq (36) is real due to Eq (35), while the last term is real by 
(k <-> k'). 

For the third order sums, we obtain in the same way 



S 3 = E ^k-ic l^o(k)| 4 [ 2|$, (k)| 2 |^ (k')| 2 - |$, (k)| 2 $: o (k) $, (k') 



k,k' 



-i^ (koi 2 ^ (k)$: o (ko 

s' 3 = 2^^, |$, (k)| 4 [ |$ n) (k')| 4 - |<Mk')| 2 *: (k)<Mk') 



(37) 
(38) 



k,k' 
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The energy shift A depending on (S3 + S'%/2) as = 0, we can check that this quantity 
is real due again to Eq (35). 

By inserting the 3D value of $> vo (k) , namely 

^ = (17^ (if • < 39 > 

into Eqs (34,36,37,38), and by calculating the integrals over k and k', we obtain 

* = f (£) E , (40) 



10795tt 3 



a 



3 



3 



= If] (42) 



with E = -e 2 /2a x . 
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FIGURE CAPTIONS 



Fig. (1) 

(a) Direct Coulomb interaction ^nij between two excitons which are scattered from 
the states to the (m, n) states, the "in" and "out" excitons being made with the 
same electron-hole pairs. 

(b) Coupling between and (m, n) excitons induced by the fact that two electrons 
and two holes can be coupled in two different ways to form two excitons. 

Fig. (2) 

Overlap A^\ as defined in Eq (12), between two iV-exciton states, when one of the 
N excitons is in a state m, instead of the ground state 0. 

(a) Diagrammatic representation of the integral equation (13) which links A^P to 
[AS^ 1 ' 1 ^ . The cross represents the A mp0 o exchange process. 

(b) Density expansion of Affl as given in Eq (15). The quantity Fn, defined in Eq 
(14), differs from 1 because the excitons are not exact bosons. 

Fig. (3) 

Overlap B^J, as defined in Eq (18), between two iV-exciton states when two of the N 
excitons are in the states (m, n) instead of the ground state 0. 

(a) Diagrammatic representation of the integral equation (20) which links Bf£2 to 
(v4.^ v ~ 2 - ) ) and (B^~ 2 ^ . The cross of the first line corresponds to A mn0 o while the two 
crosses of the last line correspond to A n<?0 o an d A mp0 o- 

(b) Density expansion of B^ as given in Eq (21). F N is defined in Eq (14). 

Fig. (4) 

(a) The sum Si defined in Eq (24) with one exchange process \ mn oo which transforms 
the (0, 0) excitons into (m, n) excitons, followed by one direct Coulomb interaction (,QQ mn 
which scatters these (m,n) excitons back into the ground state (0,0). 

(b) The sum S2, defined in Eq (25), with two exchanges and one direct Coulomb 
scattering. 

(c) (d) (e) The sums S3, S' 3 and S 3 , defined in Eq (26) and below, with three exchanges 
and one direct Coulomb scattering. 
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